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For a nonrelativistic classical particle undergoing arbitrary oscillations, the generalized effective 
potential 'I' is derived from nonlinear eigenfrequencies of the particle-field system. Specifically, the 
ponderomotive potential is extended to a nonlinear oscillator, resulting in multiple branches near 
the primary resonance. For a pair of natural frequencies in a beat resonance, scales linearly with 
the internal actions and is analogous to the dipole potential for a two-level quantum system. Thus 
cold quantum particles and highly-excited quasiclassical objects permit uniform manipulation tools, 
particularly, one-way walls. 
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I. INTRODUCTION 



II. GENERALIZED EFFECTIVE POTENTIAL 



Multiscale adiabatic dynamics of classical particles 
in oscillating and static fields is simplified within the 
oscillation-center (OC) approach, which allows separat- 
ing fast quiver motion of the particles from their slow 
translational motion [H, 0, 13| • Hence the average forces 
are embedded into the properties of the OC, yielding a 
quasiparticle with a variable effective mass Weff 0, 0] ■ In 
each given case, nies can be Taylor-expanded at nonrel- 
ativistic energies so as to appear as an effective potential 
^ e.g., ponderomotive or diamagnetic p^ . 

Yet the nonrelativistic limit must permit also an indepen- 
dent calculation of . For linear oscillations^he general- 
ized effective potential was derived in Ref. However, 
a comprehensive method of finding ^' for nonlinear quiver 
motion has not been proposed. 

The purpose of this work is to calculate, from first prin- 
ciples, the generalized effective potential ^I^ for a nonrela- 
tivistic classical particle undergoing arbitrary oscillations 
in high-frequency or static fields. We proceed by finding 
eigenmodes in the particle-field system; hence ^ is ob- 
tained like in the dressed-atom approach [111, [Hj El) [3 
but from nonlinear classical equations. Specifically, we 
show that the ponderomotive potential extended to a 
nonlinear oscillator has multiple branches near the pri- 
mary resonance. Also, for a pair of natural frequencies 
in a beat resonance, ^' scales linearly with the internal 
actions and is analogous to the dipole potential for a two- 
level quantum system. Thus cold quantum particles and 
highly-excited quasiclassical objects permit uniform ma- 
nipulation tools, particularly, stationary asymmetric bar- 
riers, or one-way walls [H [11, [13, [H, [ll, ill, HH, HI HI . 

The work is organized as follows. In SecHH we obtain 
the general form of the effective potential ^. In Sec. Illli 
we derive the equations for oscillation modes. In Sec. lIVl 
we calculate the ponderomotive potential from the in- 
finitesimal frequency shift of the oscillating field coupled 
to a particle at a primary resonance, both linear and non- 
linear; see also Appendix [X] In Sec. |Vl we find 5* near 
a beat resonance and show the analogy with the dipole 
potential. In Sec. IVIl we explain how 5" allows one-way 
walls. In Sec. IVIIl we summarize our main results. 



Consider a particle which exhibits slow dynamics in 
canonical variables (r, P) superimposed on fast oscilla- 
tions in angle-action variables ~ const), such that 
zero JT" corresponds to purely translational motion. The 
OC Lagrangian is then written as [4\ 



Co 



(1) 



where v = r, and mcfr(r, v; ^) is the effective mass, the 
dependence on being parametric [13] ■ Thus the com- 
plete Lagrangian, which describes also the oscillations, 
equals Co + 'd ■ J" , so the corresponding Hamiltonian H 
matches that of the oscillation center: 7i = P • v — £o- 

Assume nonrelativistic dynamics, i.e., v <^ c and 6m = 
rrics — m <ti m, where m is the true mass. Hence 



r 1 2 
Lo — — mv - 
2 

so the Hamiltonian reads 



U, U — Smc , 



p2 

2 m 



Zm 



(2) 



(3) 



where P = mv - d^U, and ^{J = 0) = 0. 

Following Ref. [1], one can find ^ from the relativis- 
tic particle trajectory in given fields; however, a general 
nonrelativistic approach is also possible. Formally, field 
modes can be understood as particle degrees of freedom. 
Then, like in the dressed-atom approach [ll|, [l2|, the 
effective potential can be found from eigenfrequencies 
vj = = d jTi of the particle-field system: 



■oj ■ dj. 



(4) 



In this case, ^ depends on r and P only parametrically, 
through ■cc7(r, P, J"). 

The canonical frequencies can be redefined such that 
■oj ^ vj + const, adding a constant to 5". Albeit arbi- 
trarily large, this contribution does not affect the motion 
equations, so we abandon the requirement that ^' must 



2 



remain small compared to mc^ ; hence actual physical fre- 
quencies can be used for zu. Particularly, for uncoupled 
modes one gets ^I' = 4'o, 



(5) 



where we used Eq. ^ and (fJ, J), {u>,I) for the unper- 
turbed frequencies and actions of the particle and the 
field, correspondingly. For an unbounded field (/ — > oo), 
the second term in Eq. ^ is infinite. As it is fixed 
though, the force on a particle is determined only by 
n • J and the finite modification of the effective potential 
due to coupling, 

$ = * - fo, (6) 
where $ also can be found from Eq. as shown below. 

III. PARTIAL MODE DECOMPOSITION 

Suppose weakly nonlinear oscillations both of the 
particle Q and of external fields [25] , so their Lagrangian 
reads ^) — Co + £int, where C-mt is a perturbation 
to a bilinear form Co (26l |. 

Co = lu-AH)-u-m)-li^-m- (7) 

Here M, R, Q are NxN real matrices; M and Q are sym- 
metric, R is antisymmetric, and rankM = N = dim^. 
At zero R, CqH, k) can be diagonalized to yield 



^ 1 1 



where Lj describe individual modes (27| . Then 



Qj, 



(8) 



(9) 



6 and dt standing for the variational and time deriva- 
tives. Yet, such decomposition does not hold in the gen- 
eral case, so we redefine eigenmodes, following Ref. [28 1. 
Extend the configuration space by introducing 



(e\^i-7TM-\^), |r) = (eM-V) 



(10) 



as the new, "left" and "right" , coordinate vectors, where 
TT = — R^ is the old canonical momentum. Then 



1 



Co = -[{e\m\7 



1 



{e\m\r) \ + - mr) , (11) 
where we omitted a full time derivative and introduced 

z: .XI, 5 - f ^5 *f V (12) 



M 
M 



with F = RM ^R — Q. Thus the resulting equations are 

{£\m + {e\^ = o, m\r)-^\r) = o, (13) 



both equivalent to 



Q$ = 0. 



(14) 



Eq. (dll) has 2A^ eigenmodes = e'^"^*, with i/j 
hence assumed real and nonzero; therefore, for each ^j, 
there also exists a mode — and 1^- are gener- 
ally not orthogonal. The corresponding eigenmodes of 
Eqs. ^ are 



with vector amplitudes 



and dJljk = {£j\Tl\ri;) — —2ipjk, where 



(15) 

(16) 
(17) 

(18) 



The matrix p is diagonal for distinct Vj, as seen from 
Eq. (|14p . or can be diagonalized when some of the fre- 
quencies coincide 28]; thus, 



■jk 



(19) 



Pj — —p~j. (Hence modes with Vj = are orthogonal to 
the others and can be considered separately, as implied 
below.) Therefore any {£\ and \r) are decomposed as 



(^1 



N 



N 
3 = -N 



(20) 



where the primes stand for skipping j — 0, and 



^1 = 77- {m\r,) 

zpj 



77- {e,m\r) 

zpj 



Since {£\ and \r) are real, one has rj 
tp-j = V,*; hence Eq. but with 



(21) 
■i/'j a/2 and 



The resulting equations for individual modes are 



(22) 



Pj i^j - idt) , (23) 



Particularly, at zero C 



int 7 



Jj — const. 



(24) 



Dj^Jj = (5^. /lint, 

similar to reduced Eqs. ^ 

On the other hand, Lj = [dj — Vj)Sfj] thus d^ Lj 
is also the action corresponding to the angle dj: 

J.^pMiW (25) 

so Vj is the canonical frequency. Then the mode energy 
is VjJj (thus Pj > for stable modes with Vj > 0, hence- 
forth implied), and Eq. ([5]) is recovered. 

In the next sections, we apply Eqs. ((23)l to find eigen- 
modes for nonzero Ci^t, with becoming partial os- 
cillations. Hence the effective potential modification $ 
[Eq. ©] is obtained from Eq. 
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IV. PRIMARY RESONANCE 
A. Linear oscillator 

First, we calculate for a linear coupling between a 
pair of modes ipi and -02, say, 

Ant =fTV'lV'2 +<^*V'l>2, (26) 

where a — const. In this case, Eqs. ([^5]) yields 

Dii^i = a*tP2, D2^2 = <y^i\ (27) 

hence a quadratic equation for the eigenfrequencies -ntj, 

PiP2{t^ - vi){w - V2) = \(t\^ , (28) 

from which ^' — w\J\ + VJ2J2 is obtained. 

As a particular case, consider interaction of a particle 
internal mode having frequency O and action J = p\'^^ 
with an external oscillating field E — Ee^"^^^ having fre- 
quency uj and action / = pe\E^ . Given that the field 
occupies a volume 1^ — s- 00, the frequency shifts 8^ and 
(5a; are infinitesimal, Eq. ([28| yielding 



-, 8lope = 



{Q — ll!)pe^ {uj — Q)p 

Since pE ex. V, one has dflJ <C SluI, whereas 



SluI ^ 



\E\ 



(lu — Cl)p 
is nonvanishing. Then one gets 



(29) 



(30) 



(31) 



where $0 is the so-called ponderomotive potential (for 
the general expression see Appendix [X)) . an insignificant 
constant loI is removed, and a = 4|(Tp[(0 — w)p]~^. Since 



>ii\E' 



UJ — n 



(32) 



where = |crp/p > 0, the effective potential becomes 
infinite at the linear resonance [Fig.[T]^a)]. However, non- 
linear effects remove this singularity, as we show below. 



B. Nonlinear oscillator 

Consider the effective potential near a nonlinear reso- 
nance, with a Duffing oscillator as a model system. Then 



(33) 



where f3 — const, yielding 

D^jj ^ a*E + f3\ijj\^^lj, DEE = a^jj. (34) 



Separate the driven motion from free oscillations, ip = 

-SujpEE = aX, (35) 
{n-uj)pX = a*E + 2(3\Y\^X + (3\X\^X, (36) 
-SflpY = (3\Y\^Y + 2P\X\^Y. (37) 

From Eq. it follows that SQJ - f3\XY\^, which 

we assume, for simplicity, small compared to $0 ^ SloI; 
hence VE* « flJ + <i>, where $ = /p^ Slu dl is the modi- 
fied ponderomotive potential. The field frequency shift 
is found from the cubic equation (1-1- C,h?)h = 1 for 

h = (5^//$o; C = P^o/ilP?, *o = k^lV(7P), 7 = 
UJ — n + <^(J), <; = 2/3 J/p^, and the above condition of 
negligible SflJ reads <;/(w — f2) <C 1. Thus, 



HOdc, (38) 



where /i(C) has three branches [Fig.[IJb)]. One of those is 
unstable [29]; hence two branches of the ponderomotive 
potential, $1 and <i>2- Their asymptotic form fiows from 
W^i,2(C 0±) « C and WiiC ^ 0-) » 2^ [Fig. [He)]: 



$^2 - <i>o ^ $r«2..p|i^l4M, (39) 
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and Eq. ([32]) is recovered from $^^2 ii^ the limit of small C 
and ? = [Fig.md)]. On the other hand, Wi{C ^ 00) « 
||Cp/^ so $1 ^ 3|cr£;|'*/3|/3|2/3/(2^) at the resonance, 
i.e., the singularity vanishes. 



V. BEAT RESONANCE 

Suppose that >Cint also contains a term cubic in xp, say, 
,5£int = e^iV'2*V'3 + e*V'i*V'2V'3, (40) 

where e = const. That would normally yield a potential 
small compared to $0j yet maybe except when 



1^1 ^ V2+ 1^3- 



(41) 



Below we study the beat resonance (|41]) neglecting the 
quadratic dnt, so the envelopes 'ipj — ipje^'^^*' satisfy 



-lAt 



(42) 
(43) 

-JP34 = eViiVi2*e^^*, (44) 
where A = 1^3 + 1^2 — J^i is the detuning frequency. 

A. Linear coupling 

Suppose that A is large enough, so Eqs. can 
be solved by averaging. Then split ■0^ into driven and 
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-0.1 0.1 0.2 -4t 

22/3 

f (lo-CL)IcOc 



FIG. 1: (a) Resonant ponderomotive potential $o [Eq- (|32p ] in 
units $ct = Kal-Bp/a; vs. ui — Q, in units oj. (b) Solid - squared 
normalized frequency shift h — SluI/$o vs. (, — /3$o/(7P)^, 
C,{J = 0) oc {u! — Q,)~^; dashed - approximations = 1 
and h'^ = -C^. (c) VK(C) [Eq. (|38j], consists of branches 
Wi > and W2 < 0. (d) Effective potential <I> in units 
$c = = I $(0) vs. Lo-n in units = P^^'^laEf^'^ / p, 

assuming J — and P > 0: solid - Eq. (|38|l (consists of 
branches $1 > and <I>2 < 0), dashed - Eqs. (l39|l . 



free motion: = e +Yj e "^'^j*^ where <C ?7j 



are the nonlinear frequency shifts, and 




rii = 5v2 + (5i^3 + A, 


(45) 


7/2 = (Ji^i - 5vz - A, 


(46) 


rya = 5vi - 6v2 - A. 


(47) 


Thus the corresponding equations read 






+ ^2>3), 






-V3P3X3 ^ e-YiY2 , - SvspsYs = e(YiX2 


+ XiY*), 


where insignificant oscillating terms are neg 


lected. Since 


Si^j <C A, one has 





where F = piP2P3/|ep > 0. Then, from Eqs. (|4|), one 
gets * = ^"0 + 

^ J1J2 + JiJz ^ J2J3 , .n\ 

Particularly, if the third mode corresponds to a macro- 
scopic oscillating field tuned close to the beat resonance 

cj w 17i - VI2, (49) 

then J1.2 <C / oc \E\'^ . Hence the "hybrid" ponderomo- 
tive potential that is obtained is simultaneously propor- 
tional to \E\'^ and the internal actions Jj-: 

<J> = f (J2- Ji)|^p, (50) 

where A = — (fii — $12), and kI = |ep/(piP2) > 0. 

B. Nonlinear coupling 

At A < 5vj , the internal mode equations 

-ipiiji = e*i^2Ee~'^\ -ip2^2 = ti^iE*e'^' (51) 

do not allow averaging, so the above derivation is mod- 
ified as follows. First, consider strictly periodic os- 
cillations, in which case ip = ("01, ■02)'^ rewrites as 
tp = f^C:^^yi, where fy = diag(e-^2'*/2, e'?'*/^), and 

= diag(pi, P2)'/', = arg(e^*) + tt. Then 
£ ^ . £ A 

lXl^-X2~—Xi, 1X2 = -Xi + —X2, (52) 

where e = 2k^\E\. Eqs. ([5^ yield two eigenmodes at 
frequencies ±5 A, A = (e^ + A^)^/^, so x = UT\x, where 

^^(cosQ -sine\ 

\^ sm 6 cos By ^ ' 

X = {x^,x_)'^ is constant, and 8 satisfies 

cos2e = -A/A, sin2e = e/A. (54) 

Then 

^ = fJx, (55) 

where S = C^^UC^ is a constant matrix, C-^ = 
diag(p_,_, p_)"'^/^, p± = n7± (for uniformity), 

w^ = ^{n, + n2)±^A, (56) 

X± = X±e"*''±, ^± = Tn±, and x = C'^x. 

When the oscillation parameters evolve, treat Eq. ([55]) 
as a formal change of variables for Eq. . Hence £ = 
L^. -|- L_, the field part being omitted, and 

L^ = '-f{x±xl-xlx±)-p±'^±\x±\'- (57) 
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Therefore x± ^-^e independent linear modes with frequen- 
cies G7± and conserved actions J± = p±|x±P fSec. IIIip. 
and Eq. (j4|) yields 



+ f^2)(J+ + + ! (58) 



At A > £, J± = Ji,2 for A < and J± = for A > 0, 
so Eq. ([50)1 is recovered by Taylor expansion of Eq. (|5F 



(a) 



(b) 



C. Quantum analogy 

The above classical particle is the limit of a quantum 
system with plentiful states coupled to the field simul- 
taneously and ilj being the unperturbed transition fre- 
quencies [Fig.[2la)]. Yet, with pj h, the equations (|5T|) 
are also equivalent to those describing a iwo-level system 
[3) [3 HOl) with the unperturbed eigenfrequencies Qi 
and ^2 and the Rabi frequency SIr = e [Fig.[5]^b)]. Hence 
Eq. ((58|) yields as well the dipole potential for a two-level 
quantum object, e.g., a cold atom [ll|, pp. 454-461], p^ : 



^2) 



h 
2 



A2 



(59) 



'J — IXjT are the occupation numbers (Jj 



where rij = 

hrij), satisfying + n_ = 1. Similarly, Eq. (|50l) is equiv- 
alent to the dipole potential at weak coupling ^30_, jil] : 



$ = 



1a 



(l-2ni) 



(60) 



[At ni = 0, Eq. is further analogous to Eq. ([5^ . 
with the resonance at the transition frequency ili — ri2.] 
This parallelism originates from Eq. ([23| yielding 
Schrodinger-like equations for £i„t quadratic in internal 
mode amplitudes. Therefore, for any classical particle 
governed by a hybrid potential, an analogue is possible 
which is governed by a quantum dipole potential, and 
vice versa. Hence the two types of objects permit uniform 
manipulation techniques, as also discussed in Sec. IVII 



FIG. 2: (a) Schematic of a quasiclassical particle contain- 
ing two harmonic oscillators with unperturbed frequencies Qi 
and $12- The energy spectrum for each individual oscillator 
at large energies (close-up on the right) is equidistant. The 
coupling with an oscillating field E near the beat resonance 
(|49p is described by Eqs. (ISlfl . (b) A two-level system, with 
unperturbed eigenfrequencies Qi and Q2, formally described 
by the same Eqs. (|5H) . assuming the Rabi frequency Or — e. 



Ref. [1^, such dynamics was confirmed for cyclotron- 
resonant rf fields, and, in Ref. 16], a similar scheme em- 
ploying abrupt E{z) was proposed. 

Hybrid potentials (Sec. |V| permit yet another type of 
one-way walls. Assume uniform fij and A > 0; then 
^ [Eqs. ([SS]) ] is repulsive for cold particles (Ji < 

J2) but attractive for hot particles (Ji > J2)- Thus, if 
particles incident, say, from the left are preheated (via 
nonadiabatic interaction with another field), they will be 
transmitted, whereas those cold as incident from the right 
will be repelled [Fig. [2{b)] ; hence the asymmetry. 

In agreement with the parallelism shown in Sec. IV C[ 
similar one-wa y w alls for cold atoms have been suggested 
(itI [isl. I9I 20l.l2ll [37! and enjoyed experimental verifica- 



tion [2^, |23|. However, the new result here is that quasi- 
classical particles like Rydberg atoms and molecules can, 
in principle, be manipulated in the same manner, despite 
their involved eigenspectrum is different (Fig. [2]). 



VII. CONCLUSIONS 



VI. ONE-WAY WALLS 

As an effective potential, ^ can have properties distin- 
guishing it from true potentials. Particularly, it c an y ield 



:ularly, it c an y ield 

[ilJiiJiMllii, 



asymmetric barriers, or one-way walls 11^. . 
ifl mi, HI, [li], allowing current drive [3ll33|, [Sj] and 
translational cooling pTT,, ^22,, ,3^, [11] . We explain these 
barriers as follows. 



Suppose a ponderomotive potential of the form 
[or, similarly, (|55|) ]. Given Vl — ^{z), with, say, < 0, 
the average force — —dz'^ is everywhere in +z direc- 
tion, except at the exact resonance where the effective po- 
tential does not apply [Fig. [3Ja)] . Hence particles can be 
transmitted when traveling in one direction but reflected 
otherwise, even assuming uniform E{z) [s^jlH,!!^]- In 



We propose a method to calculate the generalized ef- 
fective potential \1/ for a nonrelativistic classical parti- 
cle undergoing arbitrary oscillations in high-frequency or 
static fields. We derive ^ from the oscillation eigenfre- 
quencies in the particle-field system [Eq. (|4])], like in the 
dressed-atom approach [O, [H, [H, but from nonlinear 
classical equations [Eqs. (pS)) ]. Specifically, we show that 
the ponderomotive potential [Eqs. ([5^ . (|A6P : Fig. [TJa)] 
extended to a nonlinear oscillator has multiple branches 
near the primary resonance [Eq. (|38p : Fig. [Ijd)]. Also, 
for a pair of natural frequencies in a beat resonance, ^ 
scales linearly with the internal actions [Eqs. ((501) . fSSjl ] 
and is analogous to the dipole potential [Eqs. ((59|) . ((60)) ] 
for a two-level quantum system (Fig.jJ). Thus cold quan- 
tum particles and highly-excited quasiclassical objects 
permit uniform manipulation tools, particularly, station- 
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(a) (b) 




FIG. 3: Schematic showing two types of one-way walls (arrows 
denote particle transmission and reflection), (a) The one- 
way wall is due to a ponderomotive potential near a primary 
resonance [Eqs. (|3ip . (|32p ] with uniform d^^l < 0, so Fz = 
—dz'i^ > for all z, except at the exact resonance where the 
effective potential does not apply. [The inclined asymptote 
corresponds to = f2J + const.] (b) The one-way wall is due 
to a hybrid potential [Eqs. (|5Up . (|58p ] with uniform Qj and 
A > 0: particles incident from the left receive Ji > J2, so 
they see an attractive those incident from the right have 
Ji < J2, so they see a repulsive 'i'. 



composed of modes /i with polarizations e^. From 
Sec. IIV Al it is known that the internal energy Sft ■ J 
yields a negligible contribution to Thus 



$o = 5I^^A.^A.. (A2) 



where the infinitesimal frequency shifts Suj/^ of the field 
are found as follows. Use [3| 



A„t = 2Re[E*-d], (A3) 



where insignificant oscillating terms are removed, and d 
is the particle induced dipole moment. From Eq. (j23p . 



ary asymmetric barriers, or one-way walls (Fig. [3]). 
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APPENDIX A: GENERAL EXPRESSION FOR 
THE PONDEROMOTIVE POTENTIAL 

Consider generalization of the ponderomotive potential 
([3T|) . p2|) to multiple internal oscillations and vector field 

E = ^E^, E^ = e^^^, (Al) 



-Su^p^E^ = 4 • (A4) 

Multiply Eq. (|A4|) by i?* and substitute d = aE, where 
a is the polarizability tensor; then 



5u:^I^^^-%-a^. (A5) 



Hence summation over all field modes yields the known 
expression [ssl. l39t 



<i>o = --(E*.aE), (A6) 



which holds for any internal modes contributing to a 
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